Exit time of turbulent signals: a way to detect the intermediate dissipative range. 
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The exit time statistics of experimental turbulent data is analyzed. By looking at the exit- 
time moments (Inverse Structure Functions) it is possible to have a direct measurement of scaling 
properties of the laminar statistics. It turns out that the Inverse Structure Functions show a much 
more extended Intermediate Dissipative Range than the Structure Functions, leading to the first 
clear evidence of the universal properties of such a range of scales. 
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In stationary isotropic turbulent flows, a net flux of 
energy establishes in the inertial range, i.e. from forced 
scales, Lq, down to the dissipative scale, rj,. Energy is 
transferred through a statistically scaling-invariant pro- 
cess, characterized by a strongly non-gaussian (intermit- 
tent) activity. Understanding the statistical properties of 
intermittency is one of the most challenging open prob- 
lems in three dimensional fully developed turbulence. In 
isotropic turbulence, the most studied statistical indica- 
tors are the longitudinal structure functions, i.e. mo- 
ments of the velocity increments at distance R in the 
direction of R: 



S P (R) =< (v(x + R)-v(x)) 



- 1 p 
R > 



(1) 



Typically, one is forced to analyze one-dimensional string 
of data: the output of hot-wire anemometer. In these 
cases Taylor Frozen- Turbulence Hypothesis is used in or- 
der to bridge measurements in space with measurements 
in time. Within the Taylor Hypothesis, one has the 
large-scale typical time, To = Lq/Uq, and the dissipa- 
tive time, td = rd/Ug, where Uq is the large scale ve- 
locity field, Lo is the scale of the energy injection and 
rd is the Kolmogorov dissipative scale. As a function of 
time increment, r, structure functions assume the form: 
S p (t) =< [(v(t + t) - v(t)] p >. It is well know that 
for time increment corresponding to the inertial range, 
Td <C t <C To, structure functions develop an anomalous 
scaling behavior: S p (t) ~ t^Cp) ) where £(p) is a non linear 
function, while far inside the dissipative range, r <C Td, 
they must show the laminar scaling: S p (t) ~ t p . 
Beside the huge amount of theoretical, experimental and 
numerical studies devoted to the understanding of veloc- 
ity fluctuations in the inertial range (see [Q for a recent 
overview), only few -mainly theoretical- attempts have 
focused on the Intermediate Dissipation Range (IDR) 
By IDR we mean the range of scales, t ^ Td, 
in between the two power law ranges: the inertial and 
the dissipative range. 

The very existence of the IDR is relevant for the under- 
standing of many theoretical and practical issues. Among 
them we cite: the modelizations of small scales for opti- 
mizing Large Eddy Simulations; the possible influence of 
small scales statistics on macroscopic global quantities, 
e.g. drag-reduction due to the presence of very diluted 
polymers in the fluid B ; the validity of the Refined Kol- 
mogorov Hypothesis, i.e. the bridge between inertial and 



dissipative statistics. 

A non-trivial IDR is connected to the presence of in- 
termittent fluctuations in the inertial range. Namely, 
anomalous scaling law characterized by the exponents 
C{p), can be explained by assuming that velocity fluctua- 
tions in the inertial range are characterized by a spectrum 
of different local scaling exponents: 5 T v = v(t+T)—v(t) ~ 
T h with the probability to observe at scale r a value h 
given by P T (h) ~ T z-D{h) _ This j s the celebrated multi- 
fractal picture of the energy cascade which has been con- 
firmed by many independent experiments [EJ. The non 
trivial dissipative statistics can be explained by defin- 
ing the dissipative cut-off as the scale where the local 
Reynolds number is order of unity: 



Re{T d ) 



TdV Td 



0(1). 



(2) 



By inverting Q 
fluctuating Td'- 



we immediately obtain a prediction of a 



T d (h) 



,V(i+h) 



where for sake of simplicity we have assumed the large 
scale velocity, Uo, and the outer scale, Lq, both fixed to 
one. 

In this letter we propose, and measure in experimental 
and synthetic data, a set of new observable which are 
able to highlight the IDR properties. The main idea is 
to take a onc-dimcnsional string of turbulent data, v(t), 
and to analyze the statistical properties of the exit times 
from a set of defined velocity-thresholds. Roughly speak- 
ing a kind of Inverse Structure Functions. 
This analysis allow us to give the first clear evidence 
of non-trivial intermittent fluctuations of the dissipative 
cut-off in turbulent signals. A similar approach has al- 
ready been exploited for studying the particle separation 
statistics || . Recently, exit-time moments have also been 
studied in the time evolution of a shell model . 
The letter is organized as follows. First we define the exit- 
time probability density function and we motivate why 
this PDF is dominated by the IDR. Then, we present 
the data analysis performed in high-Reynolds number 
turbulent flows and in synthetic multi-affine signals [jlO| . 
Finally, we summarize the evidences supporting a non- 
trivial IDR and discuss possible further investigations. 
Fluctuations of viscous cut-off are particularly important 
for all those regions in the fluid where the velocity field is 
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locally smooth, i.e. the local fluctuating Reynolds num- 
ber is small. In this case, the matching between non- 
linear and viscous terms happens at scales much larger 
than the Kolmogorov scale, ~ t/~ 3 / 4 . It is natural, 
therefore, to look for observable which feel mainly lam- 
inar events. A possible choice is to measure the exit- 
time moments through a set of velocity thresholds. More 
precisely, given a reference initial time to with velocity 
v(to), we define t(Sv) as the first time necessary to have 
an absolute variation equal to Sv in the velocity data, 
i.e. \v(to) — v(to + t(Sv))\ — Sv. By scanning the whole 
time series we recover the probability density functions 
of t(Sv) at varying Sv from the typical large scale val- 
ues down to the smallest dissipative values. Positive mo- 
ments of t(Sv) are dominated by events with a smooth ve- 
locity field, i.e. laminar bursts in the turbulent cascade. 
Let us define the Inverse Structure Functions (Inverse- 
SF) as: 

E p (Sv) = <t p {Sv)> . (3) 

According to the multifractal description we suppose 
that, for velocity thresholds corresponding to inertial 
range values of the velocity differences, S Td v = v m <C 
Sv <C Vm = St v, the following dimensional relation is 
valid: 

S T v ~ r h -» t{5v) ~ Sv 1/h , (4) 

with a probability to observe a value t for the exit time 
given by inverting the multifractal probability, i.e. 

P( T „ Sv 1 " 1 ) ~ foP-DWl/* (5) 

Made this ansatz, one can write down a prediction for 
the Inverse-SF, T, p (5v) evaluated for velocity thresholds 
within the inertial range: 

Ep(fa) ~ / dh Sv^ p+3 - D ^/ h ~ 5v x{p) (6) 

where the RHS has been obtained by a saddle point es- 
timate : 

X (p)=mm{\p+3-D(h)}/h} . (7) 

h 

Let us now consider the IDR properties. 
For each p the saddle point evaluation (0) selects a partic- 
ular h = h s (p) where the minimum is reached. Let us also 
remark that from (Q) we have an estimate for the mini- 
mum value assumed by the velocity in the inertial range 
given a certain singularity h: v m (h) = 5 Td (h) v ~ v h ^ 1+h \ 
Therefore, the smallest velocity value at which the scal- 
ing Yi p (8v) ~ Sv x ^ still holds depends on both v and h. 
Namely, Sv m (p) ~ j,Mp)/i+Mp). The most important 
consequence is that for Sv < Sv m (p) the integral (||) is 
not any more dominated by the saddle point value but 
by the maximum h value still dynamically alive at that 
velocity difference, l/h(Sv) = —1 — log(V)/ log(Sv). This 
leads for Sv < Sv m (j>) to a pseudo-algebraic law: 



The presence of this p-dependent velocity range, inter- 
mediate between the inertial range, E p (Sv) ~ Sv*(p\ and 
the far dissipative scaling, Y* p (Sv) ~ Sv p , is the IDR sig- 
nature. Then, it is easy to show that Inverse-SF should 
display an enlarged IDR. Indeed, for the usual direct 
structure functions (0) the saddle point h s (p) value is 
reached for h < 1/3. This pushes the IDR to a range 
of scales very difficult to observe experimentally [|J . On 
the other hand, as regards the Inverse-SF, the saddle 
point estimate of positive moments is always reached for 
h s (p) > 1/3. This is an indication that we are probing 
the laminar part of the velocity statistics. Therefore, the 
presence of the IDR must be felt much earlier in the range 
of available velocity fluctuations. Indeed, if h s (p) > 1/3, 
the typical velocity field at which the IDR shows up is 
given by Sv m (p) ~ j/MpVU+Mp))^ that is much larger 
than the Kolmogorov value Sv rd ~ v 1 / 4 . 
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FIG. 1. Inverse Structure Functions Ei(<5u). The 
straight lines shows the dissipative range behavior (dashed) 
£i(<5v) ~ 5v, and the inertial range non intermittent behavior 
(dotted) Ei(<5«) ~ (Sv) 3 . The inset shows the direct struc- 
ture function Si (r) with superimposed the intermittent slope 
C(l) = -39. 

In Fig. 1 we plot Hi(Sv) evaluated on a string of high- 
Reynolds number experimental data as a function of the 
available range of velocity thresholds Sv. This data set 
has been measured in a wind tunnel at Re\ ~ 2000. 
Let us first make a technical remark. If one wants to 
compare the predictions (^J) and (||) with the experimen- 
tal data, it is necessary to perform the average over the 
time-statistics in a weighted way. This is due to the 
fact that by looking at the exit-time statistics we are 
not sampling the time-series uniformly, i.e. the higher 
the value of t(Sv) is, the longer it is detectable in the 
time series. Let us call ti(Sv), t%(Sv), . . . , tn(Sv) the 
string of exit time values obtained by analyzing the 
velocity string data consecutively for a given Sv. N 
is the number of times for which S T v reaches a given 
threshold. It is easy to realize |ll| that the sequential 
time average of any observable based on exit-time statis- 
tics, <t p (Sv) > t = (1 / TV) Yli=i T f 'j i s connected to the 
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uniformly-in-time multifractal average, < (•) >= 
by the relation: 



N 

<ri>(5v)>=Y,r? 

«=i 



where 7»/ 2j = i Tj takes into account the non-uniformity 
in time. Let us now go back to Fig. 1. One can see that 
the scaling is very poor. Indeed, it is not possible to ex- 
tract any quantitative prediction about the inertial range 
slope. For this reason, we have only drawn the dimen- 
sional non-intermittent slope and the dissipative slope as 
a possible qualitative references. On the other hand, (in- 
set of Fig, 1), the scaling behavior of the direct structure 
functions < |<5w(r)| >^ t^ 1 ' is quite clear in a wide range 
of scales. This is a clear evidence of IDR's contamination 
into the whole range of available velocity values for the 
Inverse-SF cases. Similar results (not shown) are found 
for higher orders S p structure functions. 
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FIG. 2. Inverse-Structure-Function Ei(<5u) versus Sv for 
the synthetic signals not smoothed (NS) and smoothed with 
time windows: ST = 4.8 • 1CT 4 , 3 • 1CT 5 , 2 ■ 1CT 6 , the straight 
line slope is obtained from the inverse multifractal prediction 

(§)• 

In order to better understand the scaling proper- 
ties of H p (Sv) we investigate a synthetic multi-affine 
field obtained by combining successive multiplications of 
Langevin dynamics [fuj . The advantage of using a syn- 
thetic field is that one can control analytically the scaling 
properties of direct structure functions in order to have 
the same scaling laws observed in experimental data. The 
signal we used is sequential in time. Therefore, it does 
not present a superposed hierarchical structure as other 
multi-affine field proposed in the past [jl2 13 1. An IDR 
can be introduced in the synthetic signals by smoothing 
the original dynamics on a moving time-window of size 
ST. Imposing a smoothing time-window is equivalent to 
fixing the Reynolds numbers, Re ~ 5T~ 4 / 3 . The purpose 
to introduce this stochastic multi-affine field is twofold. 
First we want to reach high Reynolds numbers enough 
to test the inverse- multifractal formula (R). Second, we 
want to test that the very extended IDR observed in the 
experimental data, see Fig. 1, is also observed in this 
stochastic field. This would support the claim that the 



experimental result is the evidence of an extended IDR. 
In Fig. 2 we show the Inverse-SF, Hi(5v), measured in 
the multiaffine synthetic signal at high-Reynolds num- 
bers. The observed scaling exponent, x(l), is in agree- 
ment with the prediction (Q). The same agreement also 
holds for higher moments. In Table 1, we compare the 
best fit to the Y, p (Sv) measured on the synthetic field with 
the inversion formula (Q) . As for the comparison between 
the theoretical expectation (|t]) and the synthetic data let 
us note the following points. First, in |10| it was proved 
that the signal possesses the given direct-structure func- 
tions exponents for positive moments, i.e. the £(p) expo- 
nents are in a one-to-one correspondence with the D(h) 
curve for h < 1/3. Nothing was possible to be proved 
for observables feeling the ft, > 1/3 interval and therefore 
the agreement between the inversion formula (Q) and the 
numerical results cannot be proved analytically. Second, 
because the synthetic signal is defined by using Langevin 
processes, the less singular ft-exponents expected to con- 
tribute to the saddle-point (Q) is h = 0.5. Therefore, 
the theoretical prediction, Xth(q), in Table 1 has been 
obtained by imposing h rnax = 0.5. 

Let us now go back to the most interesting question 
about the statistical properties of the IDR. In order to 
study this question we have smoothed the stochastic field, 
v(t), by performing a running-time average over a time- 
window, ST. Then we compare Inverse-SF scaling prop- 
erties at varying Reynolds numbers, i.e. for different dis- 
sipative cut-off : Re ~ ST- 4 / 3 . 

The expression (||) predicts the possibility to obtain a 
data collapse of all curves with different Reynolds num- 
bers by rescaling the Inverse-SF as follows : 

Mgg(H) HSv/u) 

bx(6T/ST ) VS ' \n{8T/8T ) ' 1 ' 

where U and ST are adjustable dimensional parameters. 
Within the same experimental (or synthetic) set up they 
are Reynolds number independent (i.e. ST independent). 
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FIG. 3. Data collapse of the Inverse-SF, T,i(8v), obtained 
by the rescaling ( |lo| ) for the smoothed synthetic signals (with 
time windows: ST = 4.8 ■ 10" 4 , 3 • 10" 5 , 2 • 10" 6 ) and the 
experimental data (EXP). The two straight lines have the 
dissipative (solid line) and the inertial range (dashed) slope. 



The rationale for the rescale (|10|) stems from the 
observation that, in the IDR, h s (p) is a function of 
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hi(5v) I ln(^) only. Therefore, identifying Re oc v" 1 , the 
relation ( |To| ) directly follows from (^). This rescaling 
was originally proposed as a possible test of IDR for di- 
rect structure functions in p| but, as already discussed 
above, for the latter observable it is very difficult to de- 
tect any IDR due to the extremely small scales involved 

Fig. 3 shows the rescaling ( pL0| ) of the Inverse-SF, T,i(5v), 
for the synthetic field at different Reynolds numbers and 
for the experimental signals. As it is possible to see, the 
data-collapse is very good for both the synthetic and ex- 
perimental signal. This is a clear evidence that the poor 
scaling range observed in Fig. 1 for the experimental sig- 
nal can be explained as the signature of the IDR. The 
same behavior holds for higher moments (not shown). 
It is interesting to remark that for a self-affine signal 
(D(h) = S(h - 1/3)), the IDR is highly reduced and the 
Inverse-SF, scaling trivially as £ p (<5v) ~ (Sv) 3p , do not 
bring any new information. 

Let us summarize the results obtained and the open 
problems. First, by defining the exit-time moments, 
Ti p (6v), we argued that they must be dominated by the 
laminar part of the energy cascade. This implies that 
they depend only on the part of D{h) which falls to the 
right of its maximum , i.e. h > 1/3. These h's val- 
ues are not testable by the direct structure functions (jl]) . 
Inverse-SF are the natural tool to test any model con- 
cerning velocity fluctuations less singular than the Kol- 
mogorov value Sv ~ r 1 / 3 . 

Second, by analyzing high-Reynolds data and synthetic 
fields, we have proved that the extension of the IDR for 
H p (5v) is magnified. Moreover, the rescaling ( |i"o| ) based 
on the assumption gives a good data collapse of all 
curves for different Reynolds numbers. This is a clear 
evidence of the IDR. 

Many questions are still open. First, the analysis of a 
wider set of experimental data could make it possible to 
quantify the agreement of the data-collapse with the pre- 
diction based on (Q) and (||). Indeed, it is easy to realize 
that, by using different parameterization for the onset of 
the viscous range, one would have predicted the existence 
of an extended IDR for 11p(6v) but with a slightly differ- 
ent rescaling procedure [p| . The quality of experimental 
data available to us is not high enough to distinguish be- 
tween the two different predictions. Analyzing different 
experimental data-sets, at different Reynolds numbers, 
could also make it possible to better explore D(h) for 
h > 1/3. This is an important question because doubts 
about the universality of these D(h) values may be raised 
on the basis of the usual energy cascade picture. For 
example, as discussed above, in the Langevin synthetic- 
data a good agreement between the multifractal predic- 
tion and the numerical data is obtained by imposing 
hmax = 0.5, similarly in true turbulent data other h rnax 
values could appear depending on the physical mecha- 
nism driving the energy transfer at large scales. 
Once the attention is focused on the exit-time statistics, 
different questions connected to the entropic properties 



of the exit-times can also be asked. For these kind of 
questions there are no a priori reasons to believe that 
the information coded in the direct-statistics is similar 
to the information coded in the inverse-statistics. Work 
is in progress in this direction. 
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TABLE I. Comparison between the Inverse-SF scaling ex- 
ponents Xsynip) measured in the synthetic signal and the in- 
version of the theoretical multifractal prediction (Q), Xthip)- 
The synthetic signal has been defined such has the D(h) func- 
tion leads to the same set of experimental f (p) exponents for 
the direct structure functions. 
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Xsyn(p) 


2.32(4) 


4.40(8) 


6.38(8) 


8.3(1) 


10.1(2) 


Xthip) 


2.32 


4.34 


6.34 


8.35 


10.35 
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